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ABSTRACT 



Context. Orion is the nearest star-forming region that hosts a significant number of young and massive stars. The energy injected by 
these OB stars is thought to have created the Eridanus superbubble. Because of its proximity, Orion is a prime target for a detailed 
investigation of the interaction between massive stars and their environment. 

Aims. We study the massive star population of Orion and its feedback in terms of energy and mass, in order to compare the current 
knowledge of massive stars with kinematic and radioactive tracers in the surrounding interstellar medium (ISM). 
Methods. We assemble a census of the most massive stars in Orion, then use stellar isochrones to estimate their masses and ages, 
and use these results to establish the stellar content of Orion's individual OB associations. From this, our new population synthesis 
code is utilized to derive the history of the emission of UV radiation and kinetic energy of the material ejected by the massive stars, 
and also follow the ejection of the long-lived radioactive isotopes -* Al and ^"Fe. In order to estimate the precision of our method, we 
compare and contrast three distinct representations of the massive stars. We compare the expected outputs with observations of -*A1 
gamma-ray signal and the extent of the Eridanus cavity. 

Results. We find an integrated kinetic energy emitted by the massive stars of 1.8;^q4 x 10^^ erg. This number is consistent with the 
energy thought to be required to create the Eridanus superbubble. We also find good agreement between our model and the observed 
^*A1 signal, estimating a mass of 5.8^2 j x 10"* Mq of -^Al in the Orion region. 

Conclusions. Our population synthesis approach is demonstrated for the Orion region to reproduce three different kinds of observable 
outputs from massive stars in a consistent manner: Kinetic energy as manifested in ISM excavation, ionization as manifested in free- 
free emission, and nucleosynthesis ejecta as manifested in radioactivity gamma-rays. The good match between our model and the 
observables does not argue for considerable modifications of mass loss. If clumping effects turn out to be strong, other processes would 
need to be identified to compensate for their impact on massive-star outputs. Our population synthesis analysis jointly treats kinematic 
output and the return of radioactive isotopes, which proves a powerful extension of the methodology that constrains feedback from 
massive stars. 

Key words. Stars: abundances, early type, winds, outflows - ISM: abundances - Gamma rays: ISM 



1. Introduction 

At a distance of only some 400 pc, the Orion region is suffi- 
• • ■ ciently close that it enables us to study its stellar population 
, 5^ ' and interstellar-gas morphology in detail. The massive star (M 
K> . > 8M0) population is dominated by the Ori on OBI association , 
5_j which includes four subgroups labelled a-d dBrown et al.ll 19941) . 
d Their ages have been estimated to range between 1 and 12 Myr. 
OB 1 is located on the near side of the den sest part of the Orion 
Molecular clouds (iMaddalena et al.lfl986h . facing the Eridanus 
cavity which extends from these molecular clouds towards the 
Sun. H a features, which coincide with a hole in the HI distri- 
bution ( Heileslll976,). toget her with X-ray emission near HI fea- 
tures dBurrows et alJfl993h outline this large interstellar cavity, 
and provide evidence for the interactions between the hot gas 
in the Er idanus cavi t y and the neutral surrounding interstellar 
medium (JBallv et al.lll991h . Many different names are used for 
th e various sub structures of Orion. We use the same convention 
as lBallvl(l2008h in his recent review. 



Feedback from massive stars plays a crucial role in the for- 
mation of stars, as it is shaping the ISM and its subsequent star 
formation activity. The main feedback originates from the ejec- 
tion of matter from massive stars through their winds and super- 
nova explosions, and from their intense emission at short wave- 
lengths into the UV. This UV radiation creates large photoion- 
ized regions around the stars, and the kinetic energy associated 
with ejection of stellar matter pushes at the IS M, together cre- 



ating large shells and cavit i es (e .g. Heiles 197^ van der Hucha 
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119871: lLeithereretaT] ll992': Maeder & Contil ll99l . Kinetic en 
ergy output and UV r adiation of massive stars were stud- 
ied in IVoss et all ( l2009h . discussing the total emission from a 
population of stars with emphasis on the differences between 
various alternative stellar models. It was found that the ki- 
netic energy from winds dominates over the supernova contri- 
bution, when integrated over the first 10 M yr after the stars 
were formed (see also iLeitherer et al. 19991). This is due to 
the high wind velocities (taken from 'Howart h & Prinial 119891 : 
Lamers et al. 1995) and high mass-loss rates o f th e most mas- 
sive stars (Castor et alJll975h and (IVink et al.ll20 00), even when 
modest wind clumping (with clumping factors of about 5) is ac- 
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counted for (iRepolust et al.ll2004t lMokiemetani2007h . To uti- 
derstand the complex interplay between the massive stars and 
their local environment, which eventually result in the evolution 
of disks in galaxies, it is necessary to first create a census of the 
radiation, energy, and matter output of individual nearby star- 
forming regions like Orion, which can be validated in terms of 
observational constraints. 

The radioactive isotope ^^Al provides an interesting inde- 
pendent view on the interaction between young stars and the 
surrounding environment. It is traced by its y-ray decay line 
at 1808.63 keV, which can be observed with y-ray telescopes. 
With a mean lifetime of ~1 Myr ^^Al is a long-term tracer of nu- 
cleosynthesis from populations of massive-star sources, as they 
eject it after synthesis in stellar cores and the supernova itself 
(iPrantzos & D iehl"l996'). Typically, massive stars eject a few 
10" ^Mq of ^^Al throu gh their winds and supernovae (SN) (see 
e.g. iLimongi & Chieffi) 2006). From y-ray observations, the to- 
talmassoff'^Al in the Milky Way is estimated to be 2.8+0.8Mq 
dPiehl et al.l 12006). Measurements of the ^^Al emission from 
Orion by the COMPTEL instrument on NAS As Compton obser- 
vatory generally confirms this scenario, with a y-ray intensity of 
~ 7.5 10"^ ph cm-2 g-i ,£,jgjjl 2002). The map of this ^''Al y-ray 
emission, though not significant in its details, shows an interest- 
ing offset of the ^^Al emission from the massive stars that are 
believed to the the source of the ^^Al, and the emission appears 
rather extended (Diehl 2002). This suggests that the radioactive 
ejecta stream into the nearby Eridanus cavity from their stellar 
association sources. A similar tracer would be the isotope ^"^Fe, 
observed in t he Galaxy globally by its 1173 keV and 1333 keV 
decay lines (Smith etal. 2005: ' Harris et al." 2005; Wang et al. 
l2007h . This isotope is presumably created in neutron capture re- 
actions in late shell burning stages of such massive stars, and 
is also emitted in the supernova explosions ('Limongi & Chie; 
120061) . ^°Fe has a mean lifetime of -3.6 Myr (Rugel et al. 200 
It has not been seen from the Orion region, which is however not 
surprising, as its y-ray intensity has been found to be ~ 15% of 
the ^^Al y-ray intensity only. 

In this paper we analyze the energy and radioactive-isotope 
output from the entire stellar content of the Orion region. We 
compare the results with observational constraints, such as the 
measured strength of the 1808.63 keV line from ^''Al decay and 
the size of the Eridanus superbubble. 



2. The massive-star content of Orion 

Star formation in the Orion region is distributed over a number 
of distinct groups. We concentrate our analysis onto the 5 most 
massive gro ups, which ar e the 4 subgroups of the Orion OB 1 
association (lBlaauwlll964l) : 

- OB la is located to the northwest of Orion's Belt region. 

- OB lb defines the belt region itself. It contains three O stars, 
^ Ori and 6 Ori, which together with the B star e Ori form 
Orion's belt, and cr Ori as another prominent member 

- OBlc partially overlaps with OBlb, extending from the Belt 
to the end of Orion's Sword. One O star, i Ori, the brightest 
star in the Sword, also belongs to this group. 

- OB Id is also called the Orion Nebular Cluster It contains two 
O stars, 6*' Ori (Trapezium) and 0^ Ori. 

- A Ori. This group is often not included in the lists of associ- 
ations with OB stars in Orion. However, it is located near the 
OB 1 association, at the Head of Orion, and at approximately 
the same distance. We therefore include it in our study. It 



contains a single O star, A Ori, after which the association is 
named. 

From published data we investigate the properties of each of 
these 5 groups individually. The results constitute assumptions 
used in our analysis, and are summarized in table [1] 

The inventory of massive stars above 2 Mq in the OB 1 as- 
sociations was analyzed in Brown et al. ( 1994), who estimated 
a total of ~610 stars in the four groups, with similar numbers 
<200 in OBla,b and c. In their table 4, they list the numbers of 
stars found and the mass ranges probed in the three groups. We 
combine these with the Salpeter mass function to estimate the 
total initial number of stars in each g roup. This gives a to tal of 
420 stars, only ~ 2/3 of the result of JBrown et alJ (Il994 . who 
used a much steeper a = 2.7 initial mass function. OB Id was 
found to host 145 stars more massive than 1 Mq (Hillenbrang 
1997), making it the smallest of the 4 OBI subgroups, an d A Ori 
contains < 50 stars above a mass of 2.5Mo CDolan & Mathieul 
1200 Ih . For compari son we convert these numbers into the 2-120 
Mq range using the lSalpeten (1195 5') mass function. 

The ages and distances of the individual clusters vary some- 
what between publications, and are n ot yet agreed on in the com- 
munit y (see e.g. the recen t reviews oflBallv'l2008l:lMuench et al] 
200i IWalteretalJl2008L JMathieu 2008;'Briceno et al. 200^ 
However, most agree that the 4 associations form a sequence 
in age and distance, with OB la being the oldest and nearest and 
OB Id the younge st and most distant (although in the study of 
iBrown et al.lll994 OBlb is significantly younger than OBlc). 
While there are large uncertainties in the absolute distances to 
the groups, the relative distances are much better understood, 
and it is therefore very unlikely that they are all at the same 
distance. This is due to the fact that the systematics affecting 
the distances are the same for the four regions. For a com- 
pilation and thoro ugh discussion of the distance studies, see 
iMuench et al.l ( 1200 8). The picture is complicated by the partial 
overlap of the groups and the possibility that they themselves 
consist of several distinct subgroups with different ages and 
distances (see e.g. 'Hardie et al 1964'; 'WaiTen & Hessei^ 119771 ; 
Guetter 1981; Gieseking 1983; Genzel & Stutzki 1989). In the 
following work we adopt a distance of ~ 410 pc to OBld as a n 
average of the three recent de terminations of . Hirota etal](l2007h ; 
Sandstrom et al. (2007); Me nten et alJ(l2007h . see e.g. the review 
of Muenchet al. (2008), which is also consistent with the re- 
sults of Jeffries (2007). OBlc is slightly closer than the OBld 
group, which places it at a distance of ~ 400 pc (Muench et al 
|2008). For OBlb we adopt a distance of ~360 pc ( Brow n et al. 
1994), and for OBla a distance of -330 pc (Briceiio et al] |2005L 
2007). We note that other relatively recent determinations find 
distances almost 100 pc further away for some of the subgroups 
(see discussions in the reviews mentioned above), and that dis- 
agreement between different methods is significantly above the 
typical ~10% errors. We adopt a distance of ~450 pc to A Ori 
dPolan & Mathieull2001[) . 

It is clear that Orio n OB 1 a i s the olde st of the groups with 
an age of -8-12 Mvr (lBlaauwll l964; W arren & Hessed 11971 
lBrownetal.lll994t iBricenoet al]i2005i) . OBld, where star for- 
mation is still underway, is cle arly the youngest g roup consist- 
ing of stars with ages 0-2 Myr dBrown et alJI 19941 ; I Hillenbrand 
1997,) . The remaining three groups have intermediate ages, but 
their exact ages are challenging to estimate. There are too few 
very high-mass stars to estimate correctly the main sequence 
turn-off mass, and the lower-mass stars have not evolved sig- 
nificantly. Age esti mates for O Blb ranges fro m 1.7+1.1 Myr 
dBrown etalFl 19941) to 8 Myr dBlaauwlll964l) . see e.g. Table 
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1 in ICaballerd (l2007l). The age of OBlc is reported to be 
in the range 3-6 Mvr (lBlaauwlll964l: IWarren & Hesseil Il977t 
[Brown et al. 1994), and A Ori is com parable to the OB lb an d 
OBlc clusters at approximately 6 Mvr dPolan & Ma thieu 200 ih. 
It is thus not clear that there is any significant age diff'erence be- 
tween these three groups. 

2. 1 . The currently most massive stars 

Previous studies of the ages of the stellar groups have focused 
mainly on less massive stars, due to their much larger numbers. 
Also the strong winds from these stars make them more chal- 
lenging to analyze, as line blanketing effects have to be taken 
into account, which has only recently become possible to do 
in detail. However, it is not clear if massive stars arrive at the 
zero-age main sequence at exactly the same time as less mas- 
sive stars. We therefore compile a list of recent determinations 
of the properties of the most massive stars in the Orion region, 
to analyze the ages and masses of the stars. The advantage of 
using these stars is that unlike lower mass stars, they move sig- 
nificantly in the log Teff - log L diagram on a timescale of Myrs. 
Given a set of evolutionary models, one can therefore derive a 
relatively precise evolutionary age even if the observational er- 
rors are large. On the other hand, the theoretical evolution of 
these massive stars is still poorly understood, and assumptions 
to derive log T^ff and log L from observations introduce relatively 
large err ors. For the Orion O B 1 associ ations, we include the stars 
listed in iBrown et all (11994 . whereas iDolan & Mathieul (1200 lb 
was used for A Ori. 

Orion hosts 7 O stars altogether. Their properties were an- 
alyzed in previous studies, using photometric data int erpreted 
using local thermodynamic equilibrium (LTE) models (iKuruczl 
1 1992 1) that were not corrected for line-blanketing. The results 
provided by such analysis are, however, very unreliable for mas- 
sive stars. We therefore derive new properties based on spec- 
troscopy rather than photometry, and using line-blanketed non- 
LTE models rather than LT E: We use the recent catalogue of 
iMaiz-Apellaniz et al.l ( |2004|) to identify the spectral typ es of the 
O star s. We then use the line-blanketed models of Marti ns et al.l 
(l2005h to estimate their effective temperatures (using the obser- 
vational scale), luminosities and surface gravities. We increased 
the sample b y adding the two v ery bright B stars (e Ori A and 
K Ori) from Searle et all (l2008h . For the Orion Nebula Cluster 
a more detailed st udy of the five most mas sive Trapezium stars 
was performed bv lSimon-Diaz et aP (l2006l) . and we use their re- 
sults for these five stars. 

In figure [T] the stars are shown in an log Teff - log L plot, 
where they are compared to four different sets of stellar evo- 
lutionary tracks and isochrones, with (upper panels and lower 
left panel) and without (lower right panel) including the ef- 
fects of rotat ion. Stellar models including rotation are taken 
from iMevnet & Maeder (2005) a nd the models without inclu - 
sion of rotation from Mevnet et al.'( 1997h: ISchaller et an(ll992h : 
iLimongi & Chieffi.(.2006.) . From this plot the masses and ages of 
the individual stars can be derived. The masses implied for the 
stars do not vary significantly between the models, but a system- 
atic shift in ages between the non-rotating and the rotating mod- 
els is apparent. For the stars in groups OBlb,c and A Ori, the av- 
erage stellar age is 0.8 Myr higher for models including rotation 
than for the ones without, whereas the five stars in group OB Id 
are on average ~ 0.7 Myr younger for the models with rotation. 
There are no significant differences between the results obtained 
from the three different non-rotating models. The results for the 
models including rotation are given in table |2l 
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Fig. 2. The evolutionary masses of the stars from table |2] versus 
their spectroscopic masses. 



For many decades O-star research has been s ubject to a se- 
vere mass discrepancy (e.g. iHerrero et al.] 11992 !). where spec- 
troscopic masses (derived from log g) and evolutionary masses 
(derived from the luminosity, and the mass-luminosity relation) 
were highly discrepant, in some cases by more than 50%. Here 
we re-investigate this issue for Orion's massive star population 
(see figure |2|. It is comforting to notice that there no longer ap- 
pears to be any significant mass-discrepancy. 

It is clear that properties of the massive stars in the groups 
OB Ib-d and A Ori are consistent with ages of a few Myr, whereas 
there are no very massive stars in OB la, in agreement with a 
higher age of ~10 Myr The 4-5 Myr ages of the stars in OBlbc 
and A Ori fall approximately in the middle of the age estimates 
derived from the less massive stars, with no evid ence for the 
much lower age of 1.7 Myr for OB lb found by Brown et al.l 
(Il994l) . The most massive star in OBld is found to be the 
youngest. The stars in OBld are all consistent with an age below 
2 Myr. 

In table [T] we summarize the assumptions on the different 
stellar groups as we use them in the following. We emphasize 
that there are considerable uncertainties on the numbers given 
in table [T] which should be evident from the discussion above 
and from figure [1] As the numbers are gathered from a large 
number of sources with varying assumptions and methodology, 
and many of these lack reliable estimates of the uncertainties, we 
have not included error estimates in the table. However, in the 
following analysis, we do estimate the sensitivity of our results 
to our assumptions. 

2.2. The population of B stars 

We also compile a list of stars with mas ses between 8 and 
20 Mq for each of the 5 subgroups (from 'Br own et al] 1 19941 : 
Dolan & Mathieu 2001; Hernandez et al. 2005), and as sign the 
assumed average age per group. We note that the use of lKuruca 
(1992) atmosphere models in these earlier studies is not quite 
appropriate for these relatively massive stars, so that the uncer- 
tainty on these estimates increases significantly. However, as the 
output of kinetic energy, matter, and UV radiation from these 
B stars is relatively small, compared to the more massive O- 
type stars, this is not important for our study. As the study of 
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Table 1. The 5 regions with massive stars in Orion. Given are the estimated number of stars above 2Mq, the age for rotating and 
non-rotating models, the highest mass possible mass for stars at these ages, the estimated number of stars that has akeady exploded 
as supernovae and the distance. 



Association Stars > 2Mo Age (rotating) Age (non-rotating) M,,,, Stars> M„p Distance 



OB la 
OB lb 
OBlc 
OB Id 
/lOri 



160 
120 
140 
60 
60 



12Myr 
5.5 Myr 
5.5 Myr 

1 Myr 
5.5 Myr 



10 Myr 


18.5 


7.3 


4.6 Myr 


45 


1.3 


4.6 Myr 


45 


1.5 


IMyr 


120 





4.6 Myr 


45 


0.65 
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450 
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Fig. 1. The 7 O stars and 2 of the bright B stars in Orion, compared to stellar main sequence tracks (solid lines) and isochrones 
(dashed lines), from four different stellar evolution models. The tracks correspond to stellar models of 12,15, 20,25,40 and 60 Mq . 
The isochrones correspond to ages between and 20 Myr with a distance of 2 Myr. The stellar models are from'Mevnet et al.l (Il997h 
(upper left),iLimongi & Chieffi (.2006.) (upper right). , Schaller et aL (.1992) (lower left) and Meynet & Maeder (2005) (lower right). 



iBrown et alj d 19941) does not give the masses of stars, we calcu- present in both the iBrown et al.l (Il994t) and [Hernandez et al.l 

(l2005h study we use the former to identify the membership of 



late these from: 

. the stellar associations, while we consider the masses given by 

M — , (1) iHernandez et al.l (l2005h to be more precise. Only two stars dif- 

AGncTshT^fi fer in masses by more than 30% between the two catalogues: 

... ^ . . ^ . HD 35439 a nd HP 37756 , which have masses of 34.3 Mq and 

where g, T,« and L are their given surface gravities, effective g g ^^ from 'Brown etalj ^^ and 11.3 M^ and 13.8 M^ in 

temperatures and luminosities per star, and G and^.b are the Hernandez et al. (2005.) . For two stars that are only present in 

gravitational and the Stetan-Boltzmann constants. For the stars 
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Table 2. The stars shown in figure [1] and their derived properties, using rotating stellar tracks. Numbers from the Henry Draper 
catalogue, commonly used names, spectral types are given, together with their spectral mass, and the initial and current masses 
assuming rotating stellar models. 



Association 



HD 



Name Spectral type Mass (Spec) Mass (current) Mass (Initial) Age 



OB lb 


36486 


S Ori A 


09.511 


17.1 


21.1 


21.4 


5.6 


OB lb 


37468 


a- Ori A 


09.5V 


15.9 


19.5 


19.6 


3.8 


OB lb 


37742 


^OriA 


09.7Ib 


20.3 


30.9 


34.3 


5.5 


OB lb 


37128 


6 Ori A 


BOIa 


45.5 


34.6 


40.8 


5.7 


OBlc 


37043 


iOriA 


09III 


19.7 


26.2 


27.3 


5.2 


OBlc 


38771 


^Ori 


B0.5Ia 


27.3 


28.0 


31.8 


6.2 


OB Id 


37022 


0' Ori C 


07Vp 


44.8 


34.7 


35.5 


1.8 


OB Id 


37041 


e^ Ori A 


09V 


27.5 


24.0 


24.1 


2.2 


OB Id 


37020 


0' Ori A 


B0.5V 


12.8 


16.1 


16.2 


2.4 


OB Id 


37023 


6/' Ori D 


B0.5V 


18.4 


17.7 


17.7 


0.6 


OB Id 


37042 


e- Ori B 


B0.5V 


9.5 


13.4 


13.4 


0.0 


^Ori 


36861 


/I Ori A 


08III 


38.8 


34.0 


37.0 


4.2 



the catalogue of lBrown et alj (1 19941) . we find unrealistically high 
masses given their spectral types. These are HD 41335 and HD 
37061 for which masses of 42.2 Mq and 50.7 Mq are found, with 
spectral types B2Vne+ and BIV. A number of the stars given in 
table |2] as well as HD 35439 mentioned above, also have unre- 
alistically high masses in that study; we ascribe this to the very 
high observational uncertainties, and assign both stars a mass of 
15 Mq in our sample. The final list of stars in the 5 regions is 
given in table |3] 

2.3. The total population of massive stars 

In addition to the observed stars, a number of massive stars are 
expected to have formerly existed in the Orion region, and ex- 
ploded as supernovae in the last 10 Myr. Knowing the age of 
an individual region and the number of stars in a given (lower) 
mass range, one can calculate the expected number of higher- 
mass stars using a distribution function for initial masses (IMF). 
In table[T]we list the expected number of exploded stars for each 
of the groups, assuming a SalpeterlMF (Salpeter 1955). 

To test consistency of our inferred stellar content of the 
groups with observations, we compare the source lists to the val- 
ues listed in table [T] We add up the number of individual stars 
above 8 Mq in table |2l table |3] and the expected number of ex- 
ploded stars from table [T] a total of 62 stars is obtained, while 
70 are expected from applying the Salpeter mass function to the 
numbers of stars above 2Mq in table [1] (excluding the 1 1 stars 
that are expected to have exploded). We note that the lower mass 
limits of the star counts in Brown et al. ( 1994) are relatively high 
(4 - 7Mq) and therefore possible errors due to deviations from 
the Salpeter law below SMp are sma ll . The extrapolation of the 
star counts from ' Hillenbrand (Il997h : iDolan & Mathieul (1200 ih 
are more uncertain as the lower mass limits in these studies were 
1 Mq and 2.5 Mq. The upper mass limit depends on the assumed 
ages of the associations and stellar evolution models. However, 
the expected number of massive stars is not sensitive to this limit, 
due to the relatively small fraction of stars at the massive end of 
the IMF. 

The only group where the number of stars in our list is sig- 
nificantly different from the IMF-expected value is A Ori, where 
4 stars above 8M0 are observed, whereas 7.8 are expected. We 
note that the mass estimates in (Dolan & Mathieu 2001) are im- 
precise for the massive stars due to inappropriate atmosphere 
models used, and that several stars are estimated just below the 
8M0 limit. Another source of bias arises from us using the num- 



ber of observed OB stars within 5 deg to estimate the richness 
of the group: a significant fraction of these could be unrelated to 
the group, inappropriately scaling up the group richness. 

As our list is compiled from various sources, applying dif- 
ferent selection criteria and analyses, it is not appropriate to use 
it to estimate the mass distribution function. Nevertheless, we 
checked if our assumed Salpeter initial mass function is com- 
patible with our stellar data: We sort the observed sources into 
4 mass bins, and compare to expectations from a Salpeter and 
a Scalo IMF. For simplicity the Scalo IMF has been normalized 
to have the same normalization at 8M0 as the Salpeter IMF, in- 
stead of deriving the normalization from the star counts of the 
individual associations. The results are shown in table |4] Only 
in the 20 - 3OM0 bin the difference is significant, and within 
uncertainties we consider both a Salpeter IMF and a Scalo IMF 
adequate to represent the overall data. 

3. Outputs from the massive stars 

We investigate the ejection of matter, of ^^Al and ^'^Fe, and the 
UV emission from the stellar groups i n Orion, us i ng the pop- 
ulation synthesis method developed bv IVoss et al.l (I2OO90 . Due 
to the proximity of Orion, the populations of stars are relatively 
well-known, as described above. We discuss three approaches 
to calculate the outputs from the stars in star-forming regions, 
comparing the results. In the first approach, the Orion stellar 
population is described by three parameters: the total number 
of stars, the average stellar age, and the age spread. The second 
(refined) approach models each of the 5 subgroups separately 
with these three parameters. In the third method, we directly use 
the observed massive stars with their parameters, together with 
estimates of those that have already exploded as supernovae. 

In our population synthesis, stellar-evolution tracks are eval- 
uated/interpolated to find the mass loss and kinetic energy from 
the stellar winds as a function of time. The supernova contribu- 
tion is added as stellar evolution terminates, assuming a canoni- 
cal ejection energy of 10^' erg. The UV emission versus time is 
found from matching stellar atmosphere models with the stel- 
lar parameters at a given time. In addition to calculating the 
cumulative stellar outputs for the entire population, we esti- 
mate the statistical deviations, caused by the random sampling 
of the IMF. Discussions of the shape of the distributions caused 
by this can be found in Cervino & Luridiana ( 2006 ); Voss et akl 
( l2009h ; lGounelle et al.l ( l2009l). We compare two ways: An ana- 
lytical formula developed by ICerviiio & Luridiana (12006!) . and 
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Table 3. The observed stars with masses > 8M0 and which are 
not listed in table |2] 



Association 



HD HIC 



Mass Mass Rer 



OB la 


35007 


25028 


8.2 


H 


OB la 


34748 


24847 


8.3 


H 


OB la 


35912 


25582 


8.5 


H 


OB la 


35762 


25493 


8.6 


H 


OB la 


36351 


25861 


9.0 


H 


OB la 


35575 


25368 


9.0 


H 


OB la 


36741 


26098 


9.3 


H 


OB la 


35777 


25480 


9.4 


H 


OB la 


35299 


25223 


9.9 


H 


OB la 


36166 


25751 


10.0 


H 


OB la 


37490 


26594 


11.1 


B 


OB la 


35439 


25302 


11.3 


H 


OB la 


35149 


25142 


11.4 


H 


OB la 


35411 


25281 


11.9 


B 


OB la 


35715 


25473 


13.1 


H 


OB la 


35039 


25044 


13.3 


H 


OB lb 


36827 


26120 


8.1 


H 


OB lb 


36779 


26106 


8.5 


H 


OB lb 


37674 


26683 


9.7 


H 


OB lb 


37744 


26713 


9.7 


H 


OB lb 


37776 


26742 


10.0 


H 


OB lb 


37903 


26816 


10.1 


H 


OB lb 


36695 


26063 


12.5 


H 


OB lb 


37479 





12.7 


B 


OB lb 


36591 


25980 


13.6 


H 


OB lb 


37756 


26736 


13.8 


H 


OBlc 


37040 


26257 


8.3 


H 


OBlc 


36629 


26000 


8.4 


B 


OBlc 


39291 


27658 


8.4 


B 


OBlc 


37209 


26345 


8.6 


B 


OBlc 


37334 


26442 


8.9 


H 


OBlc 


38051 


26908 


9.2 


H 


OBlc 


36959 


26197 


9.4 


B 


OBlc 


39777 


27929 


9.6 


H 


OBlc 


37303 


26427 


11.0 


H 


OBlc 


37018 


26237 


11.4 


B 


OBlc 


35337 


25202 


11.6 


B 


OBlc 


37481 


26535 


11.6 


B 


OBlc 


37356 


26477 


11.8 


H 


OBlc 


33328 


23972 


12.7 


B 


OBlc 


37017 


26233 


14.6 


H 


OBlc 








15.0 


B 


OBlc 


36960 


26199 


16.4 


B 


OBlc 


36512 


25923 


16.8 


B 


OBlc 


41335 


28744 


15.0 (42.2) 


B 


OB Id 


36982 





8.6 


B 


OB Id 


37061 


26258 


15.0 (50.7) 


B 


^Ori 


37232 





9.3 


D 


^Ori 


34989 





11.9 


D 


A On 


36822 





17.9 


D 














'HIPPARCOS catalogue iPerryman et al.' 19971) 

^Blgrown et al. ( 1994); H Hillenbrand 1 1997); D iDolan & Mathie j <200lh . 



Fig. 3. The Orion region modelled as a single cluster (model I) 
with a flat and a Gaussian star formation history, respectively, 
and compared to a model where the age of each subgroup is 
used (model II). From top to bottom are shown the time profiles 
of the differential and cumulative energy ejection and the amount 
of ^^Al present in the surrounding ISM. The dark and light grey 
shaded areas correspond to Icr and 2cr statistical variations of 
the flat star formation history model, determined from random 
sampling of the mass function. 



Monte Carlo sampling. Our method was f ound consistent both 
with results fr om the Starburst99 code ( Leitherer et al.ll 19991 : 
IVazquez & Leitherer 2005) . and with the results o f a similar but 
differe nt population synthesis implementation by ICerviiio et al.l 

(l2000l)^ 

In IVoss et al] (|2009|) difi'erent stellar-model inputs were an- 
alyzed and compared. In the following, we use three different 
stellar models, in order to represent the possible spread from the- 
oretical predictions. 



Table 4. The number of observed stars in 4 mass ranges, c om- 
par ed to the e xpectations according to the ISalpeted (Il955h and 
the lScald ( 1198 6) mass functions. 



Mass range 


Observed 


Salpeter 


Scalo 


>30Mo 


5 


4.8 


2.8 


20 - 30Me 


3 


7.2 


4.7 


15 - 20Mo 


8 


10.7 


8.0 


8 - 15Mo 


46 


47.7 


42.5 



genevaQS: Th e ste llar-evolution mo dels of 

iMevnet & Maeded (l2005h : iPalacios et all (l2005h including 
the effects of rotation, together with the supernova yields of 
ILimongi & Chiefli (2006). 
eneya97: The stellar-evolu tion models of iMaeder & Contil 
1994 : iMevnet et aH (Il997h without inclusion of stellar ro- 
tation, wi th enhanced mass loss , together with the supernova 
yields of I Wooslev e t al.' (1995) extracted from core sizes, 
similar to the method of Ceryino et al. ( 2000 ). 
LCQ6: The stellar-evolutio n models and supernova yields of 
iLimongi & Chieffil (l2006l) . 
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Fig. 4. The comparison of the time profiles of the kinetic energy 
ejection and the ^^Al present in the ISM for 3 different sets of 
stellar evolution models. The individual ages of the subgroups 
(model II) were used for the calculation. 



For all thr ee models we use stellar-wind velocitie s according 
to wind88 (iLamers et al.lll995tlNiedzielski & Skorzvnski 2002) 
and th e atmosMS (lKuruczll992HMartins et alJ2005HSmith et alJ 
l2005h stellar atmosphere models. 

We compare three different models for Orion's stellar popu- 
lation: 



I. Orion as one cluster: This corresponds to how we would 
have to model more distant star-forming regions, where less 
information is available on the individual groups. Given the 
estimated numbers of stars in each group from table [1] the 
whole region is expected to have 81 stars above 8 Mq. 
Results from a flat star- formation rate over the last 12 Myr 
are compared to a model in which the star-formation rate is 
Gaussian with a peak 6 Myr ago and with dispersion of 3 
Myr, where we truncate the Gaussian at 2 cr. 

II. Orion as five groups: Here we treat the 5 (4) separate groups 
individually. For each of the groups we use the parameters 
from table[T] and assume a Gaussian star-formation rate with 
a dispersion of 1 Myr, again truncated at Icr. 

III. Using the observed stars: In this model we use the observed 
stellar parameters directly, together with estimates of the 
stars that have already exploded as supernovae. We use the 
derived stellar masses directly (for the stars in table|2]we use 
the masses derived from the rotating stellar evolution mod- 
els), but assign each star the average age of the association. 
The most massive stars are generated randomly from IMF 
extrapolation, as above. 



4. Results 

When comparing models of star-forming regions to observa- 
tions, it is important to understand how our incomplete under- 
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Time [Myr] 

Fig. 5. The time profiles of the differential kinetic energy ejec- 
tion, the ^^Al and ^°Fe present in the ISM, and the emission of 
ionizing photons, from model II. The contributions from the in- 
dividual subgroups are shown. 



standing of the regions may affect such comparison. Often, the 
star formation history is poorly constrained, but may be im- 
portant considering the relatively short timescales of interest of 
~ 10 Myr In figure [3] we show the time profiles of kinetic en- 
ergy ejection from the stellar winds and supernova explosions 
and of the interstellar mass of ^^Al. Lines show the three differ- 
ent star formation histories described in section |3] flat (model 
I), Gaussian (model I) and the separate ages for each subgroup 
(model II). The results of all three models are surprisingly sim- 
ilar: Values for the current and future times are the same within 
<10%. Some differences appear in the past values, increasing 
towards the time of formation of subgroup OB Id (12 Myr ago). 
We conclude that the properties investigated in this paper are 
not sensitive to the exact star formation history for regions with 
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Fig. 6. The kinetic energy ejection and ^^ Al present in the ISM 
using model III. The figure compares the contribution from cur- 
rently observed stars with the infeiTed contributions from al- 
ready exploded stars, where the latter is divided into the wind 
contribution and the supernova contribution. 



ages above 5-6 Myr, and they cannot be used to constrain earlier 
star formation, accordingly. Shaded areas in figure [3] show the 
Icr (dark grey) and 2cr (light grey) statistical variations. These 
are derived through random sampling of the mass function (see 
e.g. Cervifio & Luridiana 2006), and are large because both the 
kinetic energy of the winds and the ejection of ^''Al strongly de- 
pend on the ZAMS mass of the stars. It is clear that these vari- 
ations are larger than the uncertainties in the star formation his- 
tory. We note that the lines indicate the average values, and that 
these probability dist ributions are stron gly asymmetric for small 
numbers of stars (see lVoss et al.ll2009l) . 

In figure |4] we compare the results of three diff'erent stellar 
evolution models and supernova yields, that are considered rep- 
resentative of the spread in theoretical predictions. For all three 
models the subgroups were modelled individually (model II). A 
spread in current values of ~ 20 - 30% can be seen, yet much 
smaller than the statistical variation. The main differences at 
early times are between the stellar evolution models including 
rotation and the ones without, with both more energy and ^^Al 
being ejected from the stars in their wind phases than from their 
supernovae. This diff'erence is mainly caused by two effects: the 
somewhat higher ages of the sub-groups infeiTed by the stellar 
models including rotation (e.g. the 2 Myr higher age of subgroup 
OB 1 a) and the enhanced wind ejection caused by stellar rotation. 

The current output is dominated by the OBlb,c and Ori A 
subgroups. At their current age (~ 4 - 6) they still contain very 
massive stars, but the most massive stars have exploded rela- 
tively recently. OB 1 a is old enough to not have any stars above 
~ 2OM0 left and is therefore mainly contributing in the current 
epoch with supernova output, whereas OB Id is so young that 



most stars have not developed strong Wolf-Rayet winds yet and 
no supernovae has exploded, so overall contributions are small. 
In figure |5] we show the time profiles of the emission of kinetic 
energy and ionizing photons, and the amount of ^^Al and ^"Fe 
present in the surrounding ISM, with the relative contributions 
from the individual subgroups. Some differences between the be- 
haviour of the different outputs can be noted. For example, the 
UV radiation from subgroup OB la has become totally insignif- 
icant, as this is linked to the most massive stars, whereas the 
energy and isotope ejection from supernovae still plays a role. In 
contrast, subgroup OB Id currently emits a high fraction of the 
total ionizing UV radiation, and some kinetic energy and ^^Al, 
but no ''°Fe, which is only ejected by supernovae. 

The results show clearly that for star-forming regions with 
< 100 massive (> SMq) stars, the random sampling of the initial 
mass function limits the physical interpretation of observations. 
For the stars still present today, the actual masses of the observed 
stars were used (model III). The solid fine in Fig. |6] represents 
their outputs and appears jagged, from actual statistical sam- 
pling. This is in contrast to the already exploded stars, where 
the average output is inferred based on the IMF. In this case 
spikes caused by supernova explosions are smeared out due to 
our lack of knowledge of the actual sampling. Clearly exploded 
stars dominate the past history of the cluster, whereas the unex- 
ploded stars dominate the future. The figure also shows that the 
current state of the system is mainly determined by the exploded 
stars. While the observations of stars are important for under- 
standing the stellar population of star-forming regions, they can 
not be used to reduce the effects of the random sampling of the 
IMF significantly. 

The choice of IMF may have some impact on the results. In 
IVoss et al] (l2009l) the contribution of the different parts of the 
mass function to the time profiles of Ekj„, UV radiation and ^^Al 
and ^"Fe were discussed (see their Fig. 13). The Scalo mass func- 
tion would decrease the number of stars above 80 Mq by a factor 
of ~ 2 and the in the 40 - 8OM0 range by ~ 40%. This will sig- 
nificantly delay and flatten the peaks of the time profiles of the 
individual groups. Our results for groups OBlb,c and A Ori will 
be affected. The results for these groups are currently dominated 
by the stars in the 30 - 5QMq range, and the values predicted by 
our models should therefore be lowered by 30-40%. OB la is old 
enough to be in the regime dominated by lower mass stars (15-30 
Mq range), where the number of stars does not differ much be- 
tween the two mass functions. The total content of massive stars 
in OB Id is known and evaluated without any IMF consideration. 

5. Comparison with observations 

We compare our model predictions against three observables 
which characterize feedback in the Orion region: the kinetic en- 
ergy as manifested in ISM excavation, the ionization resulting 
from UV output, and the radioactive materials y-ray luminosity. 
The envelopes of the massive stars are ejected through stel- 
lar winds and su pernova explosions a t typical velocities of a 
few 1000 km s~' (IWooslev et al.lll995h . and this energy creates 
large cavities around OB associations. T he flow s of supernova 
ejecta inside cavities can be very complex dMacL ow et al. 2005]), 
and the propagation might be dominated by turbulent diffusion 
from magnetic field irregularities caused by th e stellar winds 
and supernovae jParizot et al.ll2004i iBalsara & K im 2005). The 
Eridanus cavity is a typical example of a cavity in the ISM, cre- 
ated by the cumulative and sustained action of massive stars 
in the Orion OBI association. This interpretation is consistent 
with the age determinations of the OBI subgroups and the cav- 
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ity found bv lBrown et alJ (Il995h . who also estimated the energy 
required for creating the Orion-Eridanus bubble to be approxi- 
mately 1.9 X 10^^ erg. Our population synthesis model yields a 
total of 1 .8+Q 4 X 10^^ erg with ~ 40% coming from OB 1 a and the 
rest from OBib,c. OBld has not yet broken out of the surround- 
ing medium and A Ori is not connected to the Eridanus bubble. 
Our population synthesis results are therefore in agreement with 
these observations. 

The total flux of the free-free radio continuum emission can 
be used as an observable reflecting the total Lyman continuum 
luminosity of a region. As Orion covers a very large area of the 
sky (~ 600 Deg^) only a part of it has been measured by ra- 
dio telescopes. Observations of the Greater Orion Nebula (M42) 
covering the most luminous parts of Orion OB 1 d, report radio lu- 
minosities in th e 1 -25 GHz band of 300-500 Jy dPelliet al.ll993t 
Ivan der Werf & Goss 1989). This translates to a emissiori of hy - 
drogen ionizing photons of 5 - 8 x 10'*'^ ph s ' (ICondonll 19921) . 
This is significantly smaller than our estimate of ~ 3 x 10"*^ ph 
s~' on average. However, the UV output is strongly dependent on 
the most massive star in a cluster, making it highly sensitive to 
small-number statistics. Indeed the number of observed massive 
stars in OBld is significantly smaller than that from a population 
synthesis view, and thus the UV radiation is statistically very un- 
certain and almost unconstrained, with a 1-sigma confidence in- 
terval of 3x10"*^ -4x10'*'' ph s"'. The absence of stars more mas- 
sive than 45 Mq in OBld indicates that the UV radiation should 
be well below the population-synthesis predicted average found 
by integrating over the entire mass function. Indeed an integra- 
tion over the expected output from the observed stars yield an 
ionizing UV output of 10"*^ ph s"' . With a leakage of 25%-50% 
of the ionizing photons, sim ilar to what has been inferred in the 
Carina region dSmith & Bro oks 2007) this number would agree 
well with the radio continuum observations. Therefore, we do 
not consider this a significant discrepancy between predicted and 
observrd ionizing energy. 

The COMPTEL y-ray telescope has mapped the all sky dis- 
tribution of the ^^Al decay line at 1.809 MeV emission over 9 
years of observations. The results for the Orion region are pre- 
sented in Diehl (2002). Depending on the spatial model, the 
emission from the Orion region is found at a confidence level 
of 7 - 9o-, and a total flux of 2.8 - 3.7 x 10"^ ph cm^^ ^-i j^ 
found. This corresponds to a mass of ~ 4 - 5 x 10""* Mq of ^^Al 
at a distance of 400 pc, in good agreement with the results shown 
in figure[3] Calculating the emission separately for the 5 groups 
and taking into account their individual distances, we get an ex- 
pected flux of 4.5^2 ^ 10^ ph s"' from the OBI association, in 
good agreement with the observations. 

A map of the observed signal, although limited by the to- 
tal signal weakness, shows ^^Al emission in the Orion region, 
with a main peak consistent with the position of Orion OBI, and 
extended emission towards lower latitudes, suggestively aligned 
with the direction of the Orion-Eridanus bubble. Nearly all the 
flux is coming from the OBlb,c groups which are producing 
equally strong signals. A modest addition of ~ 3 x 10"^ ph s"' 
is expected to come from A Ori, which was not included in the 
observational analysis. 

6. Summary and discussion 

We analyzed the population of massive stars in the nearby star- 
forming Orion region, including the four OBI subgroups (a-d) 
and the A Ori group. We analyzed the stellar contents of the in- 
dividual groups, providing updated lists of the stars more mas- 
sive than 8 Mq. Ages of the individual groups were constrained 



based on comparison between the updated properties of the most 
massive stars and stellar isochrones. Based on these results, we 
performed a study of the ejection of kinetic energy and radioac- 
tive elements from the young massive stars in Orion. We showed 
that the current state of the region only depends modestly on the 
properties of the model, such as the star formation history and 
the stellar evolution models. Main uncertainties are due to the 
unknown population of very massive stars that exploded over 
the past 10 Myr 

The population synthesis results were compared to the en- 
ergy needed to form the Eridanus superbubble, the emission of 
hydrogen ionizing photons, and the intensity of the 1.809 MeV 
line from the decay of ^''Al, showing good agreement between 
our model estimates and the observations. The ^^Al observations 
provide a valuable tracer of the population of (now not any more 
observable) stars and thus of the cumulative action of massive 
star groups, and of the kinematics of the outflows from the mas- 
sive stars. 

Our current understanding of stellar evolution and super- 
nova models is far from complete. Different models often rely 
on similar assumptions. Showing consistency between models 
and observations is important, as it supports confidence that the 
most important effects are accounted for in models. We have 
employed different models for characterizing the Orion region's 
stellar population, and for the stellar-evolution inputs to popula- 
tion synthesis. The results show that the observed properties of 
the Orion region are consistent with these models. Differences 
among models are smaller than the statistical effects caused by 
the relatively small number o f massive star s. 

S ome recent UV studies jBouret et al.ll2003l : iFullerton et al.l 
l2006h have called for a more fundamental mass-loss rate reduc- 
tion, invoking clumping factors up to ~100, much higher than 
the currently favoured values of ~ 5, and mass-loss rate reduc- 
tions of order 10. However, other studies cast doubt on these 
conclusion s based on theoretical studies o f "macro-clumping" 
JOskinovaet al..2007 : Sundavist et alj2010l) an d emission in the 
extreme UV band (IWaldron & Cassinellil2010h . As we find good 
agreement between our population synthesis and the observa- 
tions of Orion, this could either suggest that our mass-loss rates 
are realistic (and the very large clumping factors exaggerated), 
or alternatively that some unknown process is also missing in 
the stellar models. However, we note that the wind and super- 
nova contributions to the interstellar ^''Al have not been disen- 
tangled observationally, and models with weak winds to some 
degree compensate for the low ^^ Al wind yields by having larger 
core masses and therefore p roducing higher supernova yields 
(see discussion in lLimongi & Chieffi 2006), and the production 
of ^^Al in high-clumping models have not yet been explored. 
On the other hand, we emphasize that our models are in simul- 
taneous agreement with both the kinematic and the radioactive 
tracers, which would be hard to achieve with models involving 
very large clumping factors. 

The Galaxy contains hundreds of regions of massive star- 
formation. It is important to extend our approach to other re- 
gions in order to overcome the issue of small-number statistics, 
and to further test our models. Unfortunately, many such re- 
gions are either significantly less well studied than Orion, due 
to larger distances and obscurance from the foreground, or they 
only contain modest numbers of high-mass stars. Recent stud- 
ies of the relativel y small, but nearby Scorpius-Centaurus region 
(Diehl et al.'2010') and the more distant but very massive Cygnus 
region (Martin et al. 2009) have been reported, and show overall 
agreement between the observations and models, both regard- 
ing the energetics of the regions and the ^^Al signal. A further 
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candidate target is the Carina region (ISmithI [20061) . hosting a 
large population of very young and very massive stars. Due to 
the small age, the supernova contribution to the ^''Al signal in 
this region is expected to be low. A comparison between Orion 
and Carina ^''Al signal could therefore potentially constrain the 
relative wind and supernova co ntributions, similar to what has 
been done in the Cygnus region ("Mar tin et al.ll2009l) . 

Radioactive tracers are a valuable addition to the arsenal 
of probes of star formation in the Milky Way. y-ray observa- 
tions have the potential to yield information that is complemen- 
tary to observations at other wavelengths. They are emitted on 
a timescale of Myr after star formation, and with similar decay 
timescales, they trace the cumulative action of very young stars 
in the Milky Way. As the y-rays are unaffected by extinction, the 
observations of ^^Al and ^"Fe have the potential to give a com- 
plete view of the star formation in the Milky Way, unaffected by 
the obscuring effects of the molecular clouds. Using nucleosyn- 
thesis ejecta, we can expand the studies of past activity from 
stellar groups substantially. Instrument sensitivities of current y- 
ray telescopes limit such studies to the brightest, henc e most- 
near by regions. A next generation of instruments (Grei neret alJ 
12009 ) could reach hundreds of massive-star reg ions, thus s ignifi- 
cantly extending such studies. As discussed in , Voss et al.l (l2009l) 
the correlation between ^^Al and ^"Fe can potentially eliminate 
much of the uncertainty due to small-number statistics. The life- 
times of the radioactive tracers are long enough that they can be 
carried to significant distances from the massive st ars that pro- 
duced them. The COMPTEL observations of Orion ( Diehll2002l) 
provide a hint of such displacement. On the other hand the Ufe- 
times are short enough that the radioactive elements are not uni- 
formly mixed into the ISM. They are therefore valuable tracers 
of the mixing processes in the vicinity of star forming regions. 

Acknowledgements. This research was supported by the DFG cluster of excel- 
lence 'Origin and Structure of the Universe' ^http://www.universe-cluster.dej. 



References 

Bally, J., Langer, 'W.D., Liu, ■»., 1991, ApJ, 383, 645 

Bally, J. 2008, Handbook of Star Forming Regions, Volume I: The Northern Sky 

ASP Monograph Publications, Vol. 4. Edited by Bo Reipurth, p.459, 459 
Balsara, D.S. & Kim, J., ApJ, 634, 390 

Barlow, M.J., Smith, LJ., Wilhs, A.J.,1981, MNRAS, 196, 101 
Blaauw, A. 1964, ARA&A, 2, 213 
Bouret, J.-C, Lanz, T, Hillier, D. J., Heap, S. R., Hubeny, I., Lennon, D. J., 

Smith, L. J., & Evans, C. J. 2003, ApJ, 595, 1 182 
Briceno, C, Calvet, N., Hernandez, J., Vivas, A. K., Hartmann, L., Downes, J. J., 

& Berhnd, P 2005, AJ, 129, 907 
Briceno, C, Hartmann, L., Hernandez, J., Calvet, N., Vivas, A. K., Furesz, G., 

& Szentgyorgyi, A. 2007, ApJ, 661, 1119 
Biice, C. 2008, Handbook of Star Forming Regions, Volume I: The Northern 

Sky ASP Monograph Publications, Vol. 4. Edited by Bo Reipurth, p.838, 838 
Brown, A. G. A., de Geus, E. J., & de Zeeuw, P T. 1994, A&A, 289, 101 
Brown, A. G. A., Hartmann, D., & Burton, W. B. 1995, A&A, 300, 903 
Burrows, D.N., Singh, K.P., Nousek, J.A., Garmire, G.P., Good, J., 1993, ApJ, 

406, 97 
Caballero, J. A. 2007, A&A, 466, 917 
Castor, J.I., Abbott, D.C., Klein, R.I., 1975, 195, 157 
Cervino, M., & Luridiana, V, 2006, A&A, 451, 475 
Cervino, M., Knodlseder, J., Schaerer, D., von Ballmoos, P., Meynet, G., 2000, 

A&A, 363, 970 
Condon, J. J. 1992, ARA&A, 30, 575 
Diehl, R. 2002, New Astronomy Review, 46, 547 
Diehl, R., Halloin, H., Kretschmer, K., et al., 2006, Nature, 439, 45-47 
Diehl, R., et al., in preparation. 
Dolan, C. J., & Mathieu, R. D. 2001, AJ, 121, 2124 

FeUi, M., Churchwell, E., Wilson, T. L., & Taylor, G. B. 1993, A&AS, 98, 137 
Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 
Genzel, R., & Stutzki, J. 1989, ARA&A, 27, 41 
Gieseking, F. 1983, A&A, 118, 102 



Gounelle, M., Meibom, A., Hennebelle, R, & Inutsuka, S.-i. 2009, ApJ, 694, LI 

Greiner, J., et al. 2009, Experimental Astronomy, 23, 91 

Guetter, H. H. 1981, AJ, 86, 1057 

Hardie, R. H., Heiser, A. M., & Tolbert, C. R. 1964, ApJ, 140, 1472 

Harris, M.J., Knodlseder, J., Jean, R, et al., 2005, A&A, 433, 49 

Heiles, C, 1976, ApJ, 208, L137 

Hernandez, J., Calvet, N., Hartmann, L., Briceno, C, Sicilia-Aguilar, A., & 

Berlind, P 2005, AJ, 129, 856 
HeiTcro, A., Kudritzki, R. P., Vilchez, J. M., Kunze, D., Butler, K., & Haser, S. 

1992, A&A, 261, 209 
Hillenbrand, L. A. 1997, AJ, 113, 1733 
Hirota, T, et al. 2007, PASJ, 59, 897 
Howarth, I.D., Prinja, R.K., 1989, ApJS, 69, 527 
van der Hucht.K.A., Williams, PM., The, PS., 1987, QJRAS, 28, 254 
Jeffries, R. D. 2007, MNRAS, 376, 1 109 

Kurucz, R.L., 1992, in lAU Symp. 149: The Stellar Populations of Galaxies, 225 
Lamers, H.J.G.L.M., Morton, D.C., 1976, ApJS, 32, 715 
Lamers, H.J.G.L.M., Snow, TP, Lindholm, D.M., 1995, ApJ, 455, 269 
Leitherer, C, Robert, C, Drissen, L., 1992, ApJ,401, 596 
Leitherer, C, Schaerer, D., Goldader, J.D., et al., 1999, ApJS, 123, 3 
Limongi, M., Chieffi, A., 2006, ApJ, 647, 483 

Maddalena, R.J., Morris, M., Moscowitz, J., Thaddeus, P, 1986, ApJ 303, 375 
MacLow, M., Balsara, K., Kim, J., de Avillez, M.A., 2005, ApJ, 626, 864 
Maeder, A., Conti, PS., 1994, ARA&A, 32, 227 
Maiz-Apellaniz, J., Walbom, N. R., Galue, H. A., & Wei, L. H. 2004, ApJS, 151, 

103 
Martin, P, Knodlseder, J., Diehl, R., & Meynet, G. 2009, A&A, 506, 703 
Martins, F, Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 
Mathieu, R. D. 2008, Handbook of Star Forming Regions, Volume I: The 

Northern Sky ASP Monograph Pubhcations, Vol. 4. Edited by Bo Reipurth, 

p.757, 757 
Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 
Meynet, G. & Maeder, A., 2005, A&A, 429, 581-598 
Meynet, G., Gould, M., Prantzos, N., et al., 1997, A&A, 320, 460 
Mokiem, M. R., et al. 2007, A&A, 473, 603 
Muench, A., Getman, K., Hillenbrand, L., & Preibisch, T. 2008, Handbook 

of Star Forming Regions, Volume I: The Northern Sky ASP Monograph 

Publications, Vol. 4. Edited by Bo Reipurth, p.483, 483 
Niedzielski, A., Skorzynski, W., 2002, Acta Astronomica, 52, 81 
Oskinova, L. M., Hamann, W.-R., & Feldmeier, A. 2007, A&A, 476, 1331 
Palacios, A., Meynet, G., Vuissoz, C, et al., 2005, A&A, 429, 613-624 
Panagia, N., 1973, AJ, 78, 929 
Paiizot, E., Marcowith, A., van der Swaluw, E., Bykov, A.M., Tatischeff, V, 

2004, A&A, 424, 747 
Perryman, M. A. C, et al., 1997, A&A, 323, L49 
Prantzos, N., & Diehl, R., 1996, Phys. Rep., 267, 1 
Repolust, T, Puis, J., & Herrero, A. 2004, A&A, 415, 349 
Rugel, G., et al. 2009, Physical Review Letters, 103, 072502 
Salpeter, E.E., 1955, ApJ, 121, 161 
Sandstrom, K. M., Peek, J. E. G., Bower, G. C, Bolatto, A. D., & Plambeck, 

R. L. 2007, ApJ, 667, 1161 
Scalo, J. M. 1986, Luminous Stars and Associations in Galaxies, 116, 451 
Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 
Searle, S. C, Prinja, R. K., Massa, D., & Ryans, R. 2008, A&A, 481, 777 
Simon-Diaz, S., Herrero, A., Esteban, C, & Najarro, F 2006, A&A, 448, 351 
Smith, L., Nonis, R.PF, Crowther, PA., 2005, MNRAS, 337, 1309 
Smith, N. 2006, MNRAS, 367, 763 
Smith, N., & Brooks, K. J. 2007, MNRAS, 379, 1279 
Sundqvist, J. O., Puis, J., & Feldmeier, A. 2010, A&A, 510, All 
Vacca, W.D., Garmany, CD., Shull, J.M., 1996, ApJ, 460, 914 
van der Hucht.K.A., Williams, PM., The, PS., 1987, QJRAS, 28, 254 
van der Werf, P P, & Goss, W. M. 1989, A&A, 224, 209 
Vazquez, G.A., Leitherer, C, 2005, ApJ, 621, 695 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, A&A, 362, 295 
Voss, R., Diehl, R., Hartmann, D. H., Cervino, M., Vink, J. S., Meynet, G., 

Limongi, M., & Chieffi, A. 2009, A&A, 504, 531 
Waldron, W. L., & Cassinelli, J. P 2010, ApJ, 711, L30 
Walter, F M., Sherry, W. H., Wolk, S. J., & Adams, N. R. 2008, Handbook 

of Star Forming Regions, Volume I: The Northern Sky ASP Monograph 

Publications, Vol. 4. Edited by Bo Reipurth, p.732, 732 
Wang, W., et al. 2007, A&A, 469, 1005 
Warren, W. H., Jr, & Hesser, J. E. 1977, ApJS, 34, 1 15 
Woosley, S.E., Langer, N., Weaver, T.A., 1995, ApJ, 448, 315 



10 



